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Abstract. - Taking into account stress fluctuations due to thermal noise, we study thermally 
activated irreversible crack growth in disordered media. The influence of material disorder 
on sub-critical growth of a single crack in two-dimensional brittle elastic material is described 
through the introduction of a rupture threshold distribution. We derive analytical predictions 
for crack growth velocity and material lifetime in agreement with direct numerical calculations. 
It is claimed that crack growth process is inhibited by disorder: velocity decreases and lifetime 
increases with disorder. More precisely, lifetime is shown to follow a super- Arrhenius law, with 
an effective temperature 6 — 64, where 9 is related to the thermodynamical temperature and 
9d to the disorder variance. 



Introduction. - Sub-critical rupture occurs when a material is submitted to a stress lower 
than a critical threshold. It is known since the early work of Zhurkov or Brenner [1,2] that 
sub-critical rupture can be thermally activated. However, prediction of the lifetime, i.e. the 
time it takes for a sample to break under a given load, has proved difficult when the material is 
heterogeneous. In addition, although there have been several theoretical attempts to predict 
the lifetime for homogeneous [3-5] and sometimes disordered [6, 7] elastic materials, there 
have been very few theoretical studies of the actual rupture dynamics in heterogeneous media 
when thermal activation is the driving mechanism [8-10]. Recently, this dynamics has been 
modelled for a single macroscopic crack in two-dimensional homogeneous elastic media taking 
into account stress fluctuations [11]. This model has been successfully faced with experiments 
of crack growth in fax paper sheets even though this material is actually heterogeneous [12]. 
In the present letter, we will take into account heterogeneity by introducing, in the thermal 
activation rupture model, disorder in the rupture thresholds. We will show that disorder 
slows down macroscopic crack growth contrary to previous results [8,9] that show the rupture 
dynamics to be accelerated by disorder in the case of diffuse damage. 

Model for sub-critical crack growth in homogeneous brittle media. - In ref. [11], the model 
that describes sub-critical crack growth in homogeneous brittle media assumes that thermal 
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fluctuations induce local stress fluctuations with time in the elastic material according to a 
gaussian distribution: 
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with 9 the temperature in square stress unit: 9 = YksT/V, where Y is the Young modulus, 
Ub the Boltzmann constant and T the thermodynamical temperature. In eq. I[T)). a is the 
equilibrium value of the stress a fluctuating in a volume V. We assume the volume V to 
break if the fluctuating local stress a becomes larger than a critical threshold a c . In the two- 
dimensional case where a single crack is loaded in mode 1 (uniaxial stress perpendicular to the 
crack direction), stress concentration makes almost certain that breaking occurs at the crack 
tip. Actually, at the crack tip, linear elasticity predicts a divergence of the stress. However, 
the average elastic stress a m on a volume V ~ A 3 at the crack tip does not diverge and can 
be estimated as: _ 

a m (£) = 2^2, (2) 



where £ is the crack length and o~ e the applied stress. The local stress on the crack tip volume 
V is submitted to time fluctuations according to eq. with a = a m (£) . 

The cumulative probability that the instantaneous crack tip stress is larger than the thresh- 
old o~ c is P(cr c ) = p (x) dx. Then, the average time (t w ) needed by the crack tip material 
to break on the mesoscopic scale A is directly proportional to the inverse of this cumulative 
probability: (t w ) — to/P(a c ), where to is a characteristic time for stress fluctuations. The 
statistically averaged inverse crack velocity, dt/d£, is constructed as the ratio of the average 
time (t w ) over the growth length due to rupture at each tip of the crack, dt/di = (t w )/2X. 
Then, 
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v h (a m (£),a c ,9) ' 
the crack simply propagates 



with vq = A/to- For critical crack growth, i.e. when a m (£) > cr c 
at its characteristic dynamical velocity: dtjdl — 1/vq. 

This model was successfully used to describe the growth of a single crack in a fibrous 
material such as paper [12]. In that case, the scale A corresponds to the typical paper fiber di- 
ameter. However, paper is heterogeneous and distributions in fiber size, position or orientation 
introduce some disorder which has not been taken into account in this model. 

Model extension for disordered media. - Many brittle materials are heterogeneous at a 
mesoscopic scale. Crack front roughness is one evidence for this heterogeneity. To introduce 
some disorder in the model, we assume the rupture threshold of the material, er c , to be 
distributed at the scale A with a gaussian distribution: 
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with a variance verifying v^d <C ~oZ so that the distribution can be considered normalized. For 
sake of simplicity, we will consider the crack to grow on a straight line and introduce disorder 
effects only through the threshold distribution. Then, it is clear that, during its growth, the 
crack is faced, in statistical average, with the whole intrinsic threshold distribution of the 
material. The relevance of this simplification will be discussed in the last section. 
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Using the threshold distribution of eq. Q at the crack tip, a new expression for the 
statistical average of the inverse crack velocity can be written: 



dt 
dl 



Pth(<J c )da c 
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where Vh is the crack velocity for the homogeneous case. Here, the inverse velocity, 1/Vh, is 
weighted by the threshold distribution because the waiting time for the crack to grow by a 
fixed step A is the statistical variable. Assuming Oc — a m ^> \fd + \/~9~ d , the integrand in eq. 
(0 takes significative values only for a range of <j c such that: (i) a c — a m > and the integral 
can be truncated from a m to +oo; (ii) the energy barrier (a c — a m ) 2 /29 can be considered far 
larger than one so that an asymptotic development of the complementary error function can 
be introduced in eq. JSJ. After rearrangement, we obtain: 
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Then, we can estimate the integral: 
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Here, it is important noticing that both terms of eq. @ are infinite when 8 < 9^. So, it 
is clear that, in the case where 9 < 9d, the inverse crack velocity and the rupture time will 
become infinite. On the other hand, considering the case where 9 > 9 d , we get: 
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Eq. 10 corresponds, with an additional prefactor 9/(9— 9 d ), to the asymptotic development 
of eq. for the homogeneous case introducing an effective temperature (9 c ff = 9 — 9 d . Then, 
eq. © can be seen as an approximation of: 
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when 9 > 9 d . 
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- 9 d v h (a m ,a c , 9 - 9 d ) 

The conclusion of this part is that crack velocity is lowered by the introduction of disorder 
in the model material, mainly through the appearance of a lowered effective temperature 
9 e s = 9 — 9 d in place of the thermodynamical temperature 9. The variance 9 d can then be 
seen as a temperature of disorder. 

Direct numerical analysis of the model. - In order to check the validity of eq. (|10|l . we 
study numerically the model starting from eq. (JSJ). In fig. ^l), we plot (9 — 9 d )/9 (dt/dl) nunl , 
obtained by numerical calculation of eq. (J3J)'s right hand side, as a function of the crack 
length for a set of particular values for the model parameters. The analytical law, y = 
l/vh(T m (£),o\:,9 e ff), fits very well the numerical data using 9 c h as a fitting parameter. 
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Fig. 1 - On the left: Semi-log plot of numerical (9 — 6d)/0 (dt/d£) n um as a function of £ (vo = 1, 
A = 1, Oc = 1, (T e = 0.10, 9 = 0.005, 9 d = O.OOl) and its fit by y = l/v h (a m (£), WZ, 6 aS ); on the right: 
9cg obtained from the fits as a function of 9 — 0a, with 9 and 9 d both varying (0.005 < 9 < 0.25 and 
510" 7 < d < 0.004). 



It is clear as one can see in fig. QJ)) that the effective temperature 8 e s obtained from the 
numerical data fits corresponds very well to the one predicted analytically, 9 c g — 9 — 9d, when 
9 and 9d are both varying. However, the coincidence of the numerical 9 c g with 9 — 9d is not 
perfect. Numerical 9 c h is actually perfectly described by 9 e s = a(9d) (9 — 9d) with a ranging 
from 1 to 1.014 when 9d increases from to 0.004 (ef^ = 1). These numerical results confirm 
the good quality of eq. (flCjft even if it is the result of approximate calculations. 

Rupture time. - We will now go a step further in the analysis of our model by studying 
the rupture time dependence with the model parameters. When there is no disorder, rupture 
time is denned as the time for a crack of initial length £ t to grow until it reaches a critical 
length £ c such that cr m (£ c ) = a c . In paper [11], we got an approximate expression of the mean 
rupture time for homogeneous systems assuming some approximations in eq. J2J '■ 

= t e , where Oi = a m {£i). (11) 

This law has been tested experimentally on the slow growth of a single crack in fax paper 
sheets [12]. It is shown that statistically averaged crack growth curves are in good agreement 
with the model predictions as well as the lifetime dependence on applied stress, initial crack 
length and Young modulus. 

Backing on the calculations of paper [11], we easily get an analytical expression for the 
mean lifetime in the heterogeneous case from eq. (|10f) : 

' e w- e d) = r des e ^<>-<>d) , when 9 > 9 d . (12) 



\/9 - 8 d v <j. 

We note that the lifetime is super- Arrhenius. Fig. [5] shows the logarithm of the lifetime, 
obtained by numerical integration of eq. 10 and rescaled respectively by r Q (on the left) and 
by T o es ( on the right) as a function of the energy barriers (a~ — Oi) 2 /29 (on the left) and 
(oc — a i) 2 1 '2(9 — 9d) (on the right). Eq. 111(1 does not rescale the data when 9d is varied. 
But, it is clear that the rescaling by eq. (fT2^l makes all the data for different 9 and 9d fall 
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Fig. 2 - On the left: logarithm of the rupture time as a function of the energy barrier (W^— Ui) 2 j2Q for 
many different growth conditions (0.005 < 9 < 0.25 and 5 10 -7 < Qd < 0.004); on the right: logarithm 
of the rupture time as a function of (WZ — (Ji) 2 /2(6 — 9d) for the same numerical data. 



very well on the y = x straight line. This confirms the relevance of eq. 1121 which predicts 
an increase of the lifetime with Od- Once again, analytical and numerical calculations confirm 
the deceleration of the crack growth due to disorder. 

Facing the results with previous theoretical works. - The results presented in this letter 
differ from previous theoretical works on one-dimensional Disordered Fiber Bundle Model (ld- 
DFBM) describing the influence of disorder on the fracturing process [7-10]. In this model, 
as long as the critical fraction of broken fibers is less than 50% [10], the introduction of a 
threshold distribution actually accelerates rupture with an effective temperature larger than 
the thermodynamical temperature. The rupture process starts from the weaker fibers and 
progressively affects the stronger ones. As a result, the distribution of breaking thresholds is 
truncated from below by a front moving with time towards higher thresholds. Because rupture 
can occur everywhere in the system, this dynamics is non-local and the damage rate is given 
by how many fibers break at each time step. Consequently, the effect of disorder is obtained 
by averaging the damage rate over the threshold distribution. This is not at all the process 
which occurs when a large crack propagates in a disordered system. Here, the probability of 
breaking at the crack tip is greatly enhanced by stress concentration effects even for a large 
disorder. Then, the rupture dynamics becomes local, and one needs to average over disorder 
the time to break the fiber at the crack tip, i.e. the inverse of the damage rate. Finally, it 
is worth noticing that there are similarities between our model for crack growth predicting a 
super-Arrhenius rupture time and the trap model introduced by Bouchaud [13] to describe 
dynamics in glassy systems. 

Discussion on the straight line growth hypothesis. - In our analysis, we assumed the crack 
to grow on a straight line even if roughness is experimentally observed for sub-critical cracks 
in heterogeneous materials. Is this hypothesis relevant for crack growth in heterogeneous 
materials? Two viewpoints can be raised. 

One could think that roughness is due to the fact a crack, at each moment, grow in the 
direction where the material is the weakest. In this case, the crack, during its growth, will 
not experience the intrinsic threshold distribution of the material, but only a part of this 
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distribution that corresponds mainly to the lowest thresholds. Then, the distribution to be 
used in eq. J3J is not the intrinsic threshold distribution of the material. 

But, one could also think that crack growth occurs in the direction where the stress 
intensification is the largest. Statistically, the larger the stress, the higher the probability 
to grow in the corresponding direction. This idea is supported by recent theoretical and 
experimental works [14,15]. In this framework where the crack tip stress distribution rules 
the roughness, it is fair to consider the intrinsic threshold distribution of the material to rule 
the crack growth. 

Experimental reality is probably a compromise between the two competing mechanisms: 
growth governed by stress intensification and growth following the path through the weakest 
regions. In the case of a material with a mesoscopic structure, such as paper, the crack tip is 
surrounded by a finite number of fibers so that the crack can grow only in a finite number of 
directions. The determination of the growth direction is clearly due the competition between 
stress intensification on each fiber compared to its toughness. It is worth noticing that crack 
tip structure complexity can create a complex stress distribution. Consequently, the largest 
stress is not necessarily in the crack main direction so that roughness can appear even if the 
crack always grows in the maximum stress direction. 

Some numerical work has been performed to illustrate a situation where the crack is allowed 
to "choose" between a finite number of growth directions. We model a two-dimensional elastic 
system as a network of springs forming a lattice. More details about the simulations are given 
in [11]. The only discrepancies with the simulations presented there are the use of a hexagonal 
lattice instead of the previously used square lattice and the introduction of a rupture threshold 
distribution. We can see in fig. [21 that the crack can actually "choose" between three directions 




T 



Fig. 3 - Geometry of the hexagonal spring lattice of the numerical simulation with a non-straight 
crack. 

(so called upper, straight and lower directions) at each step. In fig. 01 a typical crack pattern 
after a thermally activated growth is presented. We can notice that the crack grows almost 
straight. Actually, the stress intensification on the upper and lower springs is only about 80% 
of the one of the straight spring. So, the crack will grow through the springs on the sides only 
if the rupture threshold of the straight spring is very large. This is only a rare event and we 
can say that the crack experiences essentially the intrinsic threshold distribution truncated 
of the very large thresholds only. In this simulation, the crack path is mainly ruled by stress 
intensification. Experimentally, the number of fibers joining at the crack tip can be larger 
than three so that variations in stress intensification between fibers are probably smoother 
than in the simulation. 

No definitive conclusion on the growth process can be given because, depending on the 
considered material, one of the two competing mechanisms (growth governed by stress in- 
tensification and growth following the path through the weakest fibers) will dominate. The 
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Fig. 4 - Image of a crack created from an initial crack of length £i by simulating thermal activation 
in the hexagonal spring lattice. 

results presented in this letter hold essentially for media where the crack path is ruled by 
stress intensification. 

Conclusion. - Assuming the existence of thermal stress fluctuations in the material, we 
modelled sub-critical crack growth in disordered materials. The influence of material disorder 
on the growth of a single crack in two-dimensional brittle elastic material is described through 
the introduction of a rupture threshold distribution. Analytical predictions of the crack veloc- 
ity and material lifetime have been derived in agreement with direct numerical calculations. 
The conclusion is that the crack growth process is inhibited by disorder: velocity decreases 
and lifetime increases with the temperature of disorder 9d- The influence of disorder is simply 
accounted for by introducing an effective temperature # G ff = 9 — 9d for thermodynamical tem- 
peratures above 9d- Hence, lifetime follows a super- Arrhenius law. More experimental work is 
needed to confirm this super-Arrhenius regime. In addition, the analogy with glassy systems 
suggests that the growth of a crack in a disordered material could also occur in a "glassy" 
regime, i.e. when 9 < 9d- 

* * * 
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